On the multi-fractal structure of traded volume in financial markets 
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In this article we explore the multi-fractal properties of I minute traded volume of the equities 
which compose the Dow Jones 30. We also evaluate the weights of linear and non-linear dependences 
in the multi-fractal structure of the observable. Our results show that the multi-fractal nature of 
traded volume comes essencially from the non-Gaussian form of the probability density functions 
and from non-linear dependences. 
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I. INTRODUCTION 

The intricate character of financial markets has been one of the main motives for the physicists interest in the 
study of their statistical and dynamical properties Q, Q . Besides the asymptotic power- law behaviour for probability 
density function for price fluctuations, the return, and the long-lasting correlations in the absolute return, another 
important statistical feature observed is the return multi-fractal nature [3j . This propriety has been important in the 
establishment of analogies between price fluctuations and fluid turbulence Q and the development of multiplicative 
cascade models for return dynamics too 0. 

Changes in the price of a certain equity are basically related transactions of that equity. So, the traded volume, 
which is defined as the number of stocks that change hands during some period of time, is an important observable 
in the dynamics of financial markets. This observation is confirmed by an old proverb at Old Street that "ft takes 
volume to make prices move" @. 

In previous works several proprieties of traded volume, V , either statistical or dynamical have been studied 
In this article, we present a study of the multi-fractal structure of 1-minute traded volume time series of the 30 
equities which are used to compose the Dow Jones Industrial Average, DJ30. Our series run from the I s * of July until 
the 31 s * December of 2004 with a length of around 50fc elements each. The analysis is done using the Multi-Fractal 
Detrended Fluctuation Analysis, MF-DFA 1 Besides the multi-fractral analysis we weight the influence of 

correlation, asymptotic power-law distribution and non-linearities in the multi-fractality of traded volume. Since we 
are dealing with intra-day series we have to be cautious with the well-known daily pattern which is often considered as 
a lacklustre propriety^]. To that, we have removed that intra-day pattern of the original time series and normalised 

each element of the series by its mean value defining the normalised traded volume, v (t) — mTTjh where V (t) — SMj, 
H (r) = ' 1 N — and (. . .) is defined as the average over time (t represents the intra-day time and i the day). 

II. MULTI-FRACTALITY AND ITS COMPONENTS 

A common signature of complexity in a system is the existence of (asymptotic) scale invariance in several typical 
quantities. This scale invariance, self-affinity for time series, can be associated to a single type of structure, charac- 
terised by a single exponent, H (the Hurst exponent) or by a composition of several sub-sets, each one with a certain 
local exponent, a, and all supported onto a main structue. The former is defined as a mono-fractal and the latter as 
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1 To evaluate the local trends we have used 5 t,l -order polynomials. From this order ahead, we have a nearly polynomial-order-independent 
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FIG. 1: Left panel: Multi-fractal spectra f (a) vs. a. Right panel: Scaling exponents r(q). vs. q averaged over the 30 
companies. The legend on the right is also valid for the left. The "original" and shuffled time series present a strong multi- 
fractal character, while the shuffled plus phase randomised time series presents a narrow width in a related with the almost 
linear dependence of T B h-ran (<?) and also with the strong contribution of non-Gaussian PDF of traded volume to multi-fractality. 



a multi-fractal. In this case the statistical proprieties of the various sub-sets are characterised by the local exponents 
a are related with a fractal dimension f (<r.)|ll| composing the multi-fractal spectrum. To evaluate, numerically, this 
function we have applied the MF-DFA method [Hj. For this procedure it was proved that the g-th order fluctuation 
function, F q (s), presents scale behaviour F q (s) ~ s h ( q \ The correspondence between MF-DFA and the standard 
formalism of multifractals is obtained by, 

r(q) = qh(q)-l, (1) 

where r (q) is the exponent of the generalised partition function. From Legendre transform, / (a) = qa — r (q), we 
can relate r (q) with the Holder exponent a. Thus, using the previous equation we get 

a = h(q)+q^±, f (a) = q [a - h (q)} + 1. (2) 

In fig. njleft) we display the / (a) spectrum (full line) obtained from averages for each q over the values of the 30 
companies. In our analysis q runs from —20 to 19.5. We have verified that / (a) presents a wide range of exponents 
from a m i n = 0.32 ± 0.04 up to a max = 1.09 ± 0.04, corresponding to a deep multi-fractal behaviour. For q = 2 
we have obtained h(2) = H = 0.71 ± 0.03 which agrees with strong persistence previously observed Q. From our 
time series we can define new and related ones that can help us to quantify which factors contribute the most to the 
multi-fractal character of this observable. Among these factors we name: linear correlations, non-linear correlations 
and power-law-like PDF. To that, we have shuffled the elements (within each time series) and from these series we 
have computed h s h u f (q) ■ From these uncorrelated time series we have created another set by randomising the phase 
in the Fourier space. Afterwards, we have applied the inverse Fourier transform to come back to the time variable. 
These new series have Gaussian stationary distribution and scaling exponent h s h-ran (<?)■ 

In fig. njleft) we see that these two series also present multi-fractal spectrum, although the shuffle series has a 
wider spectrum than the shuffled plus phase randomised series. Concerning the Hurst exponent, h(2) — H, we have 
obtained H = 0.49 ± 0.03 for shuffled and H = 0.5 ± 0.03 for shuffled plus phase randomised series. Considering error 
margins, these values are compatible with H — 1/2 of Brownian motion. Furthermore, we have made a set of only 
phase randomised time series for which we have obtained H = 0.7 ± 0.03. From this values we have concluded that 
correlations have a key role in the persistent character of traded volume time series. 

Using F q (s) scaling relation for the three series ^(| an d assuming that all the factors are independent, we can 
quantify the influence of correlations, h cor (q) = h(q) — h s h u f (q), the influence of a non-Gaussian PDF, hpjyF {q) = 
h s huf (q) - h s h-ran (?), and the weight of non-linearities, h sh -ran (q) = h nUn (q) = h(q) - h cor (q) - h PDF (q). The 
multi-fractality of a time series can be analysed by mens of the difference of values between h (g max ) and h (g m in), 
hence 



Ah = h (qw) - h (g max ) 



(3) 



it is a suitable way to characterise multi-fractality. For a mono-fractal we have Aft = 0, i.e., linear dependence of 
t (g) with q. In fig. fright) we have depicted r (q) for several time series from which we have computed the various 
Aft. The results obtained are the following: Ah = 0.675, Ah cor = 0.027, Ah PDF = 0.445, and Ah rdm = 0.203. As 
it can be easily concluded the influence of linear correlations in traded volume muli-fractal nature is minimal with 
Ah cor corresponding to 4% of Aft. This value is substantially smaller than the influence of Ah n u n which corresponds 
to 30% of Aft. Our result is in perfect accordance with another previous result of us |2( where, using a non-extensive 
generalised mutual information measure |l3j , we were able to show that non- linear dependences are not only stronger 
but also more resilient than linear dependences (correlations) in traded volume time series. Last but not least, from 
the values of Aft we have verified that the main factor for the multi-fractality of traded volume time series is its 
non-Gaussian, generalised g-Gamma probability density function with a weight of 66% in Aft. Moreover, 

we have verified that the behaviour of r (q) for q > is quite different from the q < 0, which is also visible in the 
strong asymmetry of / (a) . This could indicate that large and small fluctuations appear due to different dynamical 
mechanisms. Such scenario is consistent with the super-statistical 01 approach recently presented H,@ and closely 
related with the current non-extensive framework based on Tsallis entropy [15J . Within this context and bearing in 

mind the relation 1/(1 — q sens ) = ^-/otmin — ^/ a max |l6j - we conjecture that, for traded volume, the sensivity to initial 

i 

conditions may be described by £ = [1 + (1 — q se ns) \ sens A 1 -^en.B with q sens = 0.55 ± 0.08. 
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